Dynamical effects of phonons on soliton binding in spin-Peierls systems 
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The role of dynamical magneto-elastic coupling in spin-Peierls chains is investigated by numerical 
and analytical techniques. We show that a Heisenberg spin chain coupled to dynamical optical 
phonons exhibits a transition towards a spontaneously dimerized state in a wide range of parameter 
space. The low energy excitations are characterized as solitons. No binding between solitons occurs 
in the isolated spin-phonon chain and the dynamical spin structure factor shows a broad magnon 
dispersion. However, elastic interchain coupling can lead to the formation of bound states. 
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Quasi-one-dimensional (quasi-lD) magnetic systems 
have recently received renewed experimental and the- 
oretical attention with the observation of the spin- 
Peierls transition in the two inorganic CuGeC>3 [1] and 
NaV2C>5 [2-4] compounds which consist of weakly cou- 
pled spin- 1/2 chains [5]. While such a phenomenon was 
discovered in organic materials [6], the new inorganic 
compounds can be synthesised as relatively large single 
crystals allowing for new experimental studies [1-4,7-11]. 

The phase transition was inferred from an isotropic 
drop in the magnetic susceptibility [1] at a transition 
temperature Tsp signalling a non-magnetic ground state 
(GS). The spin-Peierls transition is characterized by the 
opening of a spin gap, as has been observed by inelastic 
neutron scattering [7] (INS) and NMR spectroscopy [8,3], 
accompanied by a distortion of the lattice observed in x- 
ray diffraction experiments [9]. 

The formation at low temperature of a non magnetic 
GS has raised the possibility of observing topological 
magnetic excitations (solitons) as proposed theoretically 
in spin-1/2 frustrated Heisenberg spin chains [12] (the 
so called J\ ~ J2 model, where J\ and J2 = aJ\ are 
the nearest-neighbor and next-nearest-neighbor exchange 
couplings, respectively). At the Majumdar- Ghosh (MG) 
point, a = 0.5, the GS is doubly degenerate correspond- 
ing to the two simple dimer patterns (called A and B) ob- 
tained by a regular succession of disconnected singlet (va- 
lence) bonds. In fact, for all a > a c (a c ~ 0.241) [12-15] 
the GS was numerically shown to be dimerized. The el- 
ementary excitations in this phase are easily depicted at 
the MG point: a soliton s (antisoliton s) consists of an 
unpaired spin separating two dimer patterns A (B) and 
B (A) [16]. These objects carry spin i and can propagate 
thereby acquiring a dispersion. A spin-1 magnon excita- 
tion can be viewed as the excitation of a singlet bond into 
a triplet. However, in the J% — J2 model, such excitations 
decay into unbound soliton and antisoliton excitations. 

The spin-Peierls materials are widely described in the 
literature [17] in terms of a static ID antiferromag- 
netic frustrated dimerized Heisenberg chain [18] which 
includes, in addition to the frustrating magnetic J2 cou- 



pling, an explicit dimerization 8 of J\ [19]. Interchain 
interactions are neglected in this description. The values 
of Ji and J2 were estimated from a fit of the magnetic 
susceptibility at high temperature and the dimerization 8 
by requiring the model to have the experimental spin gap. 
Typical results such as Ji = 160 K, a = 0.36, 8 = 0.014 
were obtained for CuGe0 3 [20,14] and Ji=440 K, a = 0, 
8 = 0.048 for NaV 2 5 [4,21]. The Ji - J 2 - 8 model 
successfully predicts [22,23] the experimentally observed 
ss bound state [7,10]. Indeed, the static potential of 
strength 8 lifts the degeneracy between the two dimer 
A and B patterns. Therefore the energy cost of creating 
a A-B-A defect scales approximately as the length of 
the B structure. A confining force proportional to 8 and 
to the separation between s and s comes then naturally 
out [24-26]. Besides the spin-1 magnons, this also sug- 
gests the existence of singlet ss bound states [25,27,28] 
which could be seen e.g. in Raman spectroscopy [11]. 

The static dimerized model has, nevertheless, some 
drawbacks. It shows no spontaneous symmetry break- 
ing (since the dimerization is introduced de facto in the 
model) and ignores phonon dynamics which are expected 
to be important when the phonon frequency and the en- 
ergy scale of spin fluctuations become comparable. 

In this Letter we investigate a dynamical spin-phonon 
model which fully incorporates the phonon dynamics. 
A weak coupling renormalization group (RG) treat- 
ment [29] is used to argue that this model spontaneously 
dimerizes even in the absence of interchain coupling (at 
T = 0) for large enough spin-phonon coupling at all 
phonon frequency. Complementary numerical calcula- 
tions are needed for arbitrary parameters and show that 
this spontaneous dimerization (accompanied by a simul- 
taneous opening of a gap in the spin excitation spec- 
trum) occurs in fact in a wide range of parameter space. 
The elementary excitations are characterized as solitons. 
We show the absence of ss binding for decoupled chains. 
However, in the presence of an elastic interchain coupling 
sharp spin-1 magnon excitations are recovered in the dy- 
namical structure factor. Therefore we believe that a 
correct description of real materials must include the in- 



1 



terchain coupling [30]. 

The key ingredient of the model is the magneto-elastic 
coupling; the exchange integral is dynamically modulated 
by the relative atomic displacements along the chains. 
Using independent phonon creation (destruction) opera- 
tors b\ (6^) on each bond (i,i+l) the model reads [31,24], 

H = j£[(l + fl (& j +6|))(£.£ i+1 -±) 

i 

+ a(Si.S i+2 -^)]+H° h + H ± , 

where g is the magneto-elastic coupling constant. Here, 
we assume dispersionless optical phonons of frequency 
i.e. H° h = ft + \)- Presumably a model 

with on-site phonons will lead to similar results, how- 
ever, the model used here is probably the more relevant 
to CuGeC>3 [32]. The interchain elastic coupling H± will 
be discussed later. 

The numerical results are based on Lanczos exact diag- 
onalization (ED) [33] of closed rings with up to L=14 sites 
supplemented by finite size scaling analysis and a com- 
parison with Bcthe Ansatz exact results of the Heisen- 
berg chain and Density Matrix Renormalization Group 
(DMRG) calculations of the frustrated J\—Ji chain. A 
reliable treatment of the phonon dynamics is a difficult 
task. Preliminary studies have been carried out by con- 
sidering a single q — n phonon mode [34]. However, to 
investigate the local ss interaction it becomes necessary 
to consider the complete multi-mode problem [35]. We 
use here the variational treatment introduced by Fehren- 
bacher based on phononic coherent states. Two phononic 
states per site are retained including the vacuum and a 
phononic coherent state. This approach is non perturba- 
tive, preserves the full dynamics of the phonons, and be- 
comes exact in the weak and strong coupling limits [36] . 

As a first step, we shall consider the case of an isolated 
spin- phonon chain (H± =0). A RG argument similar 
to the one proposed in Ref. [29] can be used for this 
model for small g 2 fl/J. Due to the SU(2) symmetry 
and the absence of charge excitations there is only one 
instantaneous interaction h, with a positive bare value, 
which is marginally irrelevant [24]. Integrating out the 
phonons generates a retarded electron-electron interac- 
tion h with a bare value h(0) « —eg 2 J /ft (c is a positive 
constant) . One can then renormalize both types of inter- 
actions down to energies of the order of the phonon fre- 
quency. Below that scale, all interactions become essen- 
tially instantaneous leading to an effective Hamiltonian 
with new instantaneous interactions obtained by adding 
the retarded renormalized coupling h(£Y) to the instan- 
taneous one h(Cl). The RG equations are [29] h! = h 2 
and h! = ^hh + h 2 , involving the derivatives with respect 
to the logarithm of the energy scale. If the shift is large 
enough to change the sign of the instantaneous coupling 
(i.e. h(Q)+h(Q) < 0), the system is in the spontaneously 



dimerized phase. This scenario always occurs for a large 
enough coupling such that eg 2 > flh(0)/(2J) in the limit 

n < j. 
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figure 1: Alternation in exchange integral S* (o), spin gap 
A 01 (□) and soliton minimum energy As (0) as a func- 
tion of the coupling g for a = and frequency = 0.3 J. 

We now show numerically the existence of a T = 
spontaneous symmetry breaking towards a dimerized 
phase in a large region of parameter space. The dimer- 
ized phase is signalled, in the thermodynamic limit, by 
(i) the two fold degeneracy of the GS, (ii) the opening of 
a spin gap and (hi) a lattice dimerization. The singlet- 
triplet spin gap A 01 (L) as well as the separation A ao (L) 
between the singlet GS in the K = and K = it mo- 
mentum sectors were seen to scale accurately according 
to exponential laws [28,21]. The data strongly suggest 
A 00 = 0, a first evidence for a symmetry breaking in the 
GS. In addition, a finite extrapolated value of the spin 
gap A 01 was obtained for a wide range of parameters. As 
shown in Fig. 1, A 01 increases strongly with the magneto- 
clastic coupling g. Lower frequency phonons were found 
even more effective in opening the spin gap. 

Then, to gain more insight we have computed the 
structural distortion associated to the broken symme- 
try. The GS correlation function Ci at t(<z) = (u- q u q )o 
of the Fourier components u q = L^ 1 / 2 ^2 i exp(iqri)ui of 
the relative lattice displacements (iij = b i +b\) exhibits a 
divergence at momentum q — ir of the form CiattM oc L. 
The dimerization 5* of the exchange integral is defined 
by 5* 2 = g 2 limi^ 00 {i _1 Ci at t(7r)}. 5* obtained from a 
finite size scaling analysis is shown in Fig. 1 for a = and 
H = 0.3 J. Its g dependence follows closely the behavior 
of the spin gap. According to the previous RG analysis, 
S* should be non zero above a critical coupling, which is 
consistent with numerical results. 

We now turn to the characterization of the elementary 
magnetic excitations in the dimerized GS. If topological 
solitons exist in the spin-phonon model such excitations 
should be created in pairs. However, finite chains with an 
odd number of sites cannot accommodate a simple dimcr- 
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figure 2: (a) Minimum energy of the soliton Es as a 
function of the inverse square of the chain length 1/L 2 
for a = 0, g = (o: Bethe Ansatz), a = 0.5, g=0 (□: 
DMRG, ■: ED) and a = 0, g = 0.45, ft = 0.3J (0, ED) 
(b) Binding energy E ^ for a = 0.5, g = (□: DMRG, 
■ : ED)), a = 0, g = 0.45, ft = 0.3J (0) and a = 0.5, 
g = 0.4, ft = 0.3J (A) as a function of 1/L. 



ized pattern and the GS of such a system is expected to 
contain a single solitonic excitation. Therefore, the en- 
ergy difference defined by E s (k) = E a odd {L, k)-E° vcn (L), 
where E® dd (L, k) is the GS energy of the chain of length 
L = 2p + 1 with momentum k and E® ven (L) is the abso- 
lute GS energy of even chains interpolated at L = 2p + 1, 
can be considered as the excitation energy of a spin- 
1/2 topological defect propagating with a momentum k 
along the chain. Note that, consequently soliton ener- 
gies cannot be computed in models including an explicit 
translation symmetry breaking, like e.g. the dimerized 
Heisenberg chain. To test the accuracy of this proce- 
dure, we have compared our results for the spin-phonon 
model to the ones obtained for the purely magnetic J\- 
Ji chain at the MG point where the spin correlation 
length is the shortest. In both cases, ED of periodic 
chains show a dispersion which can be parametrized as 
E s (k) ~ \/A| + Vg(k =F tt/2) 2 , where A s is the soliton 
gap and vs a characteristic velocity. The lowest excita- 
tion in a finite system is then obtained for the momentum 
which is closest to n/2 [37]. For the spin-phonon model, 
As appears to be finite for a wide region of parame- 
ter space (see Fig. 2(a)). Hence, the spin-1/2 excitation 
spectrum is massive indicating the existence of solitons 
contrary to the case of the Heisenberg chain where spin- 
1/2 excitations (called spinons) are gaplcss [38], as explic- 
itly shown in Fig. 2(a) using the Bethe Ansatz solution 
for a = [39]. The soliton gap should not depend on 
boundary conditions, and using DMRG data for open 
chains [37], we find A s ~ 0.1170(2)J at the MG point in 
the absence of a spin-phonon coupling in good agreement 
with the estimate in Ref. [16] of As — 0.125 J. 

Spin-1 magnons can be considered as a ss combination 



in a triplet state. Consequently, the excitation triplet 
and second singlet gaps are written as A ob ~ 2Ag + E^, 
where i?ff accounts for a finite range spin-dependent in- 
teraction between s and s. Ef/? is finite only if the ss com- 
bination forms a bound state, otherwise it vanishes and s 
and s separate to infinity. The latter scenario occurs for 
example in the spontaneously (a > a c ) dimerized frus- 
trated Heisenberg chain [24] where A 00 = A 01 . However, 
in the case of spin-phonon models, a direct comparison 
between A 00 and A 01 requires some caution: low energy 
spin-0 excitations of phononic character are likely to exist 
and are indistinguishable from the magnetic singlet ss ex- 
citations. The triplet binding energy E^ = A 01 — 2Ag 
can nevertheless be calculated on closed rings and ED 
results are shown in Fig. 2(b). The data for the spin- 
phonon model are very similar to the DMRG data on 
open chains at the MG point (also shown in Fig. 2(b)) 
strongly suggesting a vanishing binding energy for all pa- 
rameters. This is consistent with the weak coupling RG 
treatment: one obtains an effective field theory similar 
to the one describing the pure spin system (with renor- 
malized coupling constants) and no bound state is ex- 
pected [24]. We conclude that solitons and antisolitons 
are not bound in the ID spin-phonon model. 




figure 3: Dynamical spin factor structure S zz (q,to) cal- 
culated at momenta q = 2im/L, n = 0, ■ • ■ , L/2 (L=12) 
(a) for a = 0, g = 0.3, ft = 0.2 J and no interchain cou- 
pling; (b) for a = 0, g = 0.08, ft = 0.2J and A = 0.1J. 
High energy peaks are due to finite size effects. A small 
broadening e = 0.04J was used. 

Because of the decay of the ss pair, we do not ex- 
pect any sharp J-function structure in the low frequency 
spin structure factor. The latter which can be exper- 
imentally obtained by INS is defined as S zz (q,Lu) = 
£„ \{<S> n \S z {q)\$ )\ 2 b-{LO -E n + E Q ), where E (E n ) is 
(are) the energy(ies) of the GS <I>o (triplet states $„) and 
S z (q) is the Fourier transform of S l z . Results for the spin- 
phonon model on a 12 site chain are shown in Fig. 3(a) 
for realistic parameters of NaV 2 5 (i.e. leading to a spin 
gap A 01 ~ 0.2 J) and are strikingly different from those 
obtained for a static dimerization [21]. The main struc- 
ture, although reminiscent of the well-known spinon dis- 
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persion of the Heisenberg chain (u>(q) — ^J\s'mq\) [38], 
is shifted towards higher energies and is much broader. 
Furthermore, the relative weight of the low energy peak 
(e.g. at q = n/2) decreases for increasing system size. 
Consistent with the previous analysis, there are no well 
defined spin-f magnons. 

Next, we argue that the interchain coupling is crucial 
in order to produce well defined magnon excitations as 
experimentally observed. To test this scenario we have 
considered an elastic interchain coupling of the form [24] , 

Hj_ = K± J2i X)<7, 7 '> u l u l > where (7,7') refers to adja- 
cent chains. We shall treat this interchain elastic cou- 
pling at the mean field level [31,24], while retaining the 
full dynamics of the ID phonons. A given chain 7 will 
then be clastically coupled to the static deformation 
(u] }o = (— l)*Mo of the Z neighboring chains. Conse- 
quently we get an additional dynamical term, -Hj_ ; mf = 
^Ejt -1 )'^; + &I) where A = ZK±u which explicitly 
doubles the unit cell. 

Turning on a small mean field interchain coupling A in 
the absence of frustration could again favor one of the 
two lattice distortions and confine solitons into pairs, in 
which case a number of ss stable bound states propor- 
tional to 1/ A (in the limit A <C J) would be expected [24] . 
Although the previous analysis on soliton binding can- 
not be done anymore, we give further arguments in favor 
of soliton binding based on the study of the dynamical 
structure factor. This one is shown on Fig. 3(b) for pa- 
rameters relevant to NaV2C>5 with a non-zero value of 
A. Qualitatively, the low energy magnon structure is 
now much better defined and the magnon dispersion is 
clearly apparent. The maximum of the dispersion occurs 
atwsi 1.7 J, an energy close to the exact value f J of 
the Heisenberg chain [38]. At momentum q = ir/2 where 
finite size effects are shown to be quite small [23,21], the 
average frequency (ui) and the width Auj of the struc- 
ture depend strongly on the presence of an interchain 
coupling. The structure clearly gets narrower when the 
interchain coupling increases. Furthermore, when A ^ 
the relative weight of the lowest energy peak increases 
with the length of the chain, contrary to the case of the 
isolated chain. 
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